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Abstract 

Interference between the photonic dipole operator Q 7 and the current-current operators Qig 
gives one of the most important QCD corrections to the B —>• X s 7 decay rate. So far, the O (a 2 ) 
part of this correction has been known in the heavy charm quark limit only (m c mj/ 2 ). 
Here, we evaluate this part at m c = 0, and use both limits in an updated phenomenological 
study. Our prediction for the CP- and isospin-averaged branching ratio in the Standard Model 
reads = (3.36 ± 0.23) x 10 -4 for E 1 > 1.6 GeV. 


1 Introduction 


The inclusive weak radiative decay 5 —> X s 7 is known to provide valuable tests of the Standard 
Model (SM), as well as constraints on beyond-SM physics. Measurements of its CP- and isospin- 
averaged branching ratio B S1 at the T(4S') experiments, namely CLEO [1], Belle M and 
Babar USKZj, contribute to the following world averagfQ [ 8 ] 

B \:* p = (3.43 ± 0.21 ± 0.07) x 10 “ 4 (1.1) 


for 5 7 > Eo = 1.6 GeV in the 5-meson rest frame. A significant suppression of the experimental 
error is expected once Belle II begins collecting data in a few years from now [10,11]. 

Let us describe the relation of B S1 to decay rates in an untagged measurement at T(4S'). 
One begins with the CP-averaged decay rates 

p T(B° -> AVy) + T(B° -> AVy) T(B~ -> AVy) + T(B* -» AVy) 

Their isospin average Y = (r 0 ± r±)/2 and asymmetry A 0 ± = (To — r±)/(r 0 + T±) are 
related to B S1 as follows 


B 


t b oT 


(1 + ryry 
\ l+r f 


± A 0 ± 


1 — rfr T \ 
l + r f ) ' 


(1.3) 


Here, ry = t b +/t b o = 1.076 ±0.004 [ 8 ] and r/ = f + ~/f 00 = 1.059 ±0.027 [ 8 ] are the measured 
lifetime and production rate ratios of the charged and neutral 5-mesons at T(4S'). The term 
proportional to A 0 ± in Eq. (11.31) contributes only at a permillc level, which follows from the 
measured value of A 0 ± = —0.01 ± 0.06 (for 5 7 > 1.9 GeV) [71 IT21[T3] . 

The final state strangeness in Eq. (11.21) (—1 for X s and ±1 for A g ) as well as the neutral 5- 
rneson flavours have been specified upon ignoring effects of the 5°5° and K°K° mixing. Taking 
the K°K° mixing into account amounts to replacing X s and Xg by A| s | with an unspecified 
strangeness sign, which leaves To and T-t invariant. Next, taking the 5°5° mixing into account 
amounts to using in T 0 the time-integrated decay rates of mesons whose flavour is fixed at the 
production time. Such a change leaves r 0 practically unaffected because mass eigenstates in 
the 5°5° system are very close to being orthogonal ( \p/q\ = 1) and having the same decay 
width [13]. In the following, we shall thus ignore the neutral meson mixing effects. 

Theoretical calculations of the 5 —> A" s y decay rate are based on the equality 


P(5 —* X s ^) Ei>Eq — r(6 —> Xg'y) E ^ >Eo ± ^Tnonp, (1-4) 

where the first term on the r.li.s. stands for the perturbatively calculable inclusive decay rate of 
the b quark into charmless partons Xf = s, sg, sgg, sqq,... and the photon. For appropriately 
chosen E 0 , the second term 5r non p becomes small, and is called a non-perturbative correction. 
For E 0 = 1.6 GeV, the uncertainty due to poor knowledge of 5r non p has been estimated to 

1 The new semi-inclusive measurement by Belle [9j which supersedes [2] is not yet taken into account in this 
average. 
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remain below 5% of the decay rate [2j. The non-perturbative correction is partly correlated 
with the isospin asymmetry because <5r nonp depends on whether B = B° or B = B~ [15]. 

As far as the perturbative contribution T(6 —>■ is concerned, its determination with an 

accuracy significantly better than 5% is what the ongoing calculations aim at. For this pur¬ 
pose, order O(a^) corrections need to be evaluated. Moreover, resummation of logarithmically 
enhanced terms like (a s ln(M^/m^))" is necessary at each order of the usual a s -expansiono 
Such a resummation is most conveniently performed in the framework of an effective theory 
that arises after decoupling of the electroweak-scale degrees of freedom. In the SM, which we 
restrict to in the present paper, one decouples the top quark, the Higgs boson and the gauge 
bosons W ± and Z°. Barring higher-order electroweak corrections, all the relevant interactions 
are then described by the following effective Lagrangian: 

8 2 

v;,v a V cmq, + V QWiQi - Q’t) . 

i =1 i=l 

(1.5) 

where Gp is the Fermi constant, and V vj are the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements. The operators Q^ are given by 

Qi = (s Llll T a u L )(u L ^T a b L ), 

Ql = (s L .l lt u L ){u L 'fb L ), 

Q i = (s L ^T a c L )(c L rT a b L ), 

Q 2 = (s L 'y fl c L )(c L ^b L ), 

Q3 = (s L 'y ll b L )^(q^q), 
q 

Qi = (sLl^bQYi^l^q), 

q 

Q 5 = {s L -f fll 'y fl2 -f fl3 b L )^2(q^ 1 Y 2 'y fl3 q), 

q 

Qe = 

q 

Qi = 

Os = (1.6) 

where the sums in Q 3,...,6 g° over ah the active flavours q = u,d, s, c, b in the effective theory. 

Decoupling (matching) calculations give us values of the electroweak-scale Wilson coefficients 
Ci(fio), where /i 0 (.). Next, renormalization group equations are used to evolve them 

down to the low-energy scale, i.e. to find Ci(nb), where /x& ~ mb/2 is of order of the final 

2 After the resummation, subsequent 0(1), 0(a s ) and 0(a 2 ) terms in this expansion are called Leading 
Order (LO), Next-to-Leading Order (NLO) and Next-to-Next-to-Leading Order (NNLO). 


Geff — ^QCDxQED^, d, S, C, b) + 


4 G f 

~V2 
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hadronic state energy in the 5-meson rest frame. Determination of the Wilson coefficients 
Ci,...,8(/h>) up to O(ag) in the SM was completed in 2006 [T51 - H9] . Matching calculations up to 
three loops [T6] and anomalous dimension matrices up to four loops [19] were necessary for this 
purpose. The three-loop matching calculation has recently been extended to the Two-Higgs- 
Doublet-Model case m■ Most of the final results have been presented for the so-called effective 
coefficients 


cfUi) = 


Ci(fjL), for 2 = 1,. 

C 7 (n) + Ej=i for i = 7 , 

C%{n) + Ej =1 ZjCjin), for i = 8, 


.., 6 , 


(1.7) 


where the numbers y 3 and Zj are such that the LO decay amplitudes for b —> sy and b —>■ sg are 
proportional to the LO terms in CfHjM,) and C| ff (fib), respectively [21] . In the MS scheme with 
fully anticommuting y 5 , one finds y = (0, 0, — — |, —y, —y ) and z = (0, 0, 1, — 20, — y ) [22] . 

Once the Wilson coefficients have been found up to the NNLO, one proceeds to 

evaluating all the on-shell decay amplitudes that matter at this order foiHl 


r(i-t^)«,>a, = G ^°“ iy,;W E cfM x 


i,j =1 


X 


G ( °\e 0 ) + + (^) 2 G|f(So.w) + 0(al) 


+ ...,( 1 . 8 ) 


where ellipses stand for higher-order electroweak corrections. At the LO, the symmetric matrix 
Gf? takes the form 


Gf(E„) = S„S J7 + 7f, 


(1.9) 


where T-p describe small tree-level contributions to b —> sqq'y from Q“ 2 and Qz,...fi p3l.;2^]. 
At the NLO and NNLO, numerically dominant effects come from Gp -), Gpp and Gpp. While 
Gj 7 is known in a complete manner [2771129] . calculations of G\ 7 and G y 27 are still in progress. 
Contributions from massless and massive fermion loops on the gluon lines have been found 
in Refs. [3QH32] . and served as a basis for applying the Brodsky-Lepage-Mackenzie (BLM) 
approximation [33]. The remaining (non-BLM) parts of G^ 2 ^ 7 have been known so far in the 
heavy charm quark limit only ( m c 'G> m^j 2) [3TU35J . 

In the present work, we evaluate the full GG C 7 for m c = Eq = 0. It is achieved by calculating 
imaginary parts of several hundreds four-loop propagator-type diagrams with massive internal 
lines. Next, both limits are used to interpolate in m c those parts of the non-BLM contributions 
to Gp\^ 7 whose exact m c -dependence is not yet known. It will give us an estimate of their 
values at the measured value of m c , and for non-vanishing E 0 . 

3 Following the notation of Ref. [25], we use tilde over G in the r.h.s. of Eq. to indicate the overall 
normalization to ml le . 
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Our current approach differs in several aspects from the one in Ref. [31] where interpolation 
in m c was applied to a combined non-BLM effect from all the with i,j £ {1,2,7, 8 }@ In 
the present paper, the only interpolated quantities are the above-mentioned parts of 2 ^ 7 . 

Exact m c -dependence of most of the other important non-BLM contributions to G\y is now 
available thanks to calculations performed in Refs. [29 [ l32 [ l36], Last but not least, the current 
analysis includes the previously unknown m c -independent part of G\% [37], all the relevant 
BLM corrections to G\f with i, j y 7 [3T ] l38 [ l39]. tree-level contributions T® [231124] . four- 
body NLO corrections m, as well as the updated non-perturbative corrections mmm- 
The only contributions to G^ with i , j £ {1,2, 7, 8 } that remain neglected are the unknown 
(n > 3)-body final state contributions to the non-BLM parts of G^ with i,j y 7. 

The article is organized as follows. In Section [2] we describe the calculation of G^ 12 j 7 for 
m c = Eq = 0. A new phenomenological analysis begins in Section [3] where m c -dependence 
of the considered correction is discussed, and the corresponding uncertainty is estimated. In 
Section [4] we evaluate our current prediction for B S1 in the SM, which constitutes an update 
of the one given in Ref. m We conclude in Section [5] Appendix A contains results for all the 
massless master integrals that were necessary for the calculation in Section [2] Several relations 
to quantities encountered in Ref. [43] are presented in Appendix B. In Appendix C, we collect 
some of the relevant NLO quantities. Appendix D contains a list of input parameters for our 
numerical analysis together with a correlation matrix for a subset of them. 


2 Calculation of G$ and for m c = £o = 0 

2.1 The bare calculation 


Typical diagrams that had to be evaluated for the present project are shown in Fig. [0 They 
represent a subset of possible unitarity cut contributions to the 6 -quark self-energy due to 
the interference of various effective operators. At the highest loop level, i.e. four-loops, this 
interference involves the operators Q \,2 and Q 7 . We need to consider two-, three- and four- 
particle cuts. Possible five-particle cuts would necessarily involve real cc pairs originating 
from the Q \,2 operator vertices, while open charm production is not included in B —>■ AT ,7 
by definition. For this reason, we skip the diagrams with five-particle cuts together with all 
the diagrams with real cc production or virtual charm loops on the gluon lines. In Section [3] 
contributions from virtual charm loops on the gluon lines will be taken over from the m c 7 ^ 0 
calculation of Ref. [32], and added to the final result. 

For efficiency reasons, we work directly with cut diagrams and employ the technique first 
proposed in [44]. The idea of the method is to represent cut propagators as 


—27T iS(p 2 — m 2 ) 


1 1 

p 2 — m 2 + ie p 2 — m 2 — ie 


( 2 . 1 ) 


4 At the NNLO level, we neglect the small Wilson coefficients C 3 ,... ,Cq, and the CKM-suppressed effects 
from Q'i 2 - 
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Figure 1: Sample diagrams for G 


( 2 ) 

(1,2)7 


with some of the possible cuts indicated by the dashed lines. 


As long as we perform only algebraic transformations on the integrands, there is no difference 
between the first and second terms on the r.li.s. of the above equation, and it is sufficient 
to work with one of them only. This is particularly convenient for the integration-by-parts 
(IBP) method for reduction of integrals [55J. The only difference in such an approach between 
complete integrals and cut integrals is that a given integral vanishes if the cut propagator 
disappears due to cancellation of numerators with denominators. This fact reduces the number 
of occurring integrals in comparison to a computation without cuts. 

In practice, the calculation follows the standard procedure. Diagrams are generated with 
DiaGen [56], the Dirac algebra is performed with FORM [57], and the resulting scalar integrals 
are reduced using IBP identities with IdSolver [56] . The main challenge of this calculation 
begins after these steps. The amplitudes for the interference contributions are expressed in 
terms of a number of master integrals, most of them containing massive internal 6-quark lines 
and a non-trivial phase space integration in D — 4 — 2e spacetime dimensions, with up to four 
particles in the final state. A feeling for the size of the problem can be gained from Tab. [D 

Having a large number of massive cut integrals, it is advantageous to devise a strategy to treat 
them in a uniform manner. It is clear that purely massless cut integrals are easier to calculate 
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n D 

nos 

^eff 

^massless 

two-particle cuts 

292 

92 

143 

9 

three-particle cuts 

267 

54 

110 

11 

four-particle cuts 

292 

17 

37 

7 

total 

851 

163 

290 

27 


Table 1: Number of diagrams no, number of massive on-shell master integrals nos, number of 
effectively computed massive master integrals n e s, and number of massless master integrals n mass i ess . 
The last two columns are explained in the text. 

than massive ones. Therefore, we aim at replacing a calculation of massive propagator integrals 
by a calculation of massless ones. This can be achieved by extending the integral definitions. 
We assume, namely, that the external momentum squared pi is a free parameter, and treat 
coefficients X) in the e-expansion of the master integrals as functions of a single dimensionless 
variable x = pi/ml- IBP identities give us differential equations 

-^(Y) = , (2.2) 

3 

with Jij(x) being certain rational functions of x. Boundary conditions for these equations in 
the vicinity of x = 0 are given by asymptotic large-mass expansions, i.e. by power-log series 
in x. A few leading terms in the series for each X t can be found by calculating products of 
massive tadpole integrals up to three loops and massless propagator ones up to four loops, as 
illustrated in Fig. [2j Next, higher-order terms can be determined from the differential equations 
themselves by substituting 2 ] in terms of power-log series in x. For our application it turns out 
that around 50 terms are sufficient to obtain the desired accuracy. This gives us high-precision 
boundary conditions at small but non-vanishing x for solving the differential equations 02 .21) 
numerically. 

On the way from the vicinity of x = 0 to the physical point at x = 1, one often encounters 
spurious singularities on the real axis. To bypass them, the differential equations are solved 
along ellipses in the complex x plane. Several such ellipses are usually considered to test whether 
the numerical solution is stable. 

Naively, one might think that as long as there are no infinities at x = 1, the numerical 
solution could be continued up to that point. However, there is an essential singularity there, 
and the integrals behave as (1 — x) n ln m (l —x), with n, m > 0 being some positive powers. Due 
to such a behaviour, the numerical solution has poor convergence, as the algorithms assume 
locally polynomial behaviour of the considered functions. In order to overcome this problem, we 
perform another power-log expansion around x — 1, and match it onto the numerical result. To 
determine the maximal power of the logarithms, we begin with observing that the highest poles 
in the cut diagrams could potentially be of order 1/e 6 , due to the presence of collincar and soft 
divergences. The coefficient of the leading singularity contains no ln(l — x) because logarithms 
are generated by expanding expressions of the form (1 — x) ae /e 6 (with a being some constant) 
in the framework of expansion by regions. Thus, finite parts of the master integral expansions 
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Figure 2: Diagrammatic representation of the asymptotic large mass expansion of two non-planar 
master integrals.Thick and thin lines represent massive and massless propagators, respectively, while 
dashed lines show the unitarity cuts. 


may only contain ln 6 (l — x). Higher powers may be needed due to the presence of spurious 
singularities, i.e. poles in the coefficients at the master integrals in the physical amplitude. In 
practice, we have used an ansatz with logarithm powers up to fifteen. Our numerical matching 
has shown that such high powers never occur in the considered problem, i.e. the respective 
expansion coefficients are consistent with zero to very high numerical precision. Using the 
matched series, we finally obtain the required values of the original master integrals at x = 1 . 
The solution procedure is schematically represented in Fig. [3^. 

Since the master integrals are considered for x ^ 1, their overall number n e s is larger than 
it would be for x — 1, i.e. n e ff > nos■ However, the massless integrals that are necessary to 
determine the boundary conditions near x = 0 are not only simpler, but also their number 
Umassiess is much smaller than nos, as seen in Tab. [0 All the massless integrals that we had to 
consider are depicted in Appendix A, in Fig. [7] and Tab. [3j 

Using the above method, we have obtained the following bare NNLO results for the consid¬ 
ered interferences in the Feynman-’t Hooft gauge: 

P^(2)bare _ 

Lr 1? — 

7^(2)bare _ 

^27 — 

+ 


1 p;f2)bare 80 1592 + 54 tt 2 

- q G ?7 + g^2 +- 2437 - + 42 - 0026519628 > 


4 30332 + 4327r 2 

3e3 2187 e 2 


67.66077706444119 


nni 


32 

729 e 


+ 0.6520676315J + m 


352 


+ 


+ 44.5070537274 
11624 228656 188 


729e 2 2187e 


+ 


6561 


243 


-7T 
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Figure 3: Left (a): Integration contour in the complex x plane. The numerical integration (Nl) is 
performed between the regions close to x = 0 and x — 1 that are accessible by power-log expansions 
(PLE). Right (b): Diagrams that give the terms marked with k in Eq. (12.31) . 


n b 


352 
729 e 


+ 


5.17409838118169 


+ 15.1790288135 I + 0(e). 


(2.3) 


Here, ni and n b denote numbers of massless and massive (m = m b ) quark flavours, while k — 1 
marks contributions from the diagrams in Fig. [3 ]d describing interferences involving four-body 
sqq'y final states and no ccy couplings. The terms proportional to rq and n b but not marked by 
k reproduce (after renormalization) the m c —> 0 limits of what is already known for non-zero 
m c [30H32] , For compactness, all the results in this subsection are given for /x 2 = e 7 m^/( 4n), 
where 7 is the Euler-Mascheroni constant. 

Some of the numbers in Eq. (12.31) have been given in an exact form even though our calcula¬ 
tion of the master integrals at x = 1 is purely numerical. However, the accuracy is very high (to 
around 14 decimals), so identification of simple rationals is possible. Moreover, renormalization 
gives us relations to lower-order results where more terms are known in an exact manner (see 
below). For the Tp-term, after verifying numerical agreement with Refs. [301139] . we have made 
use of the available exact expressions]))] Several other numbers in this subsection that have been 
retained in a decimal form can actually be related to quantities encountered in Ref. [43], as 
described in Appendix B. 

Let us now list all the lower-order bare contributions that are needed for renormalization. 
For this purpose, it is convenient to express Eq. (1 1. 8 j) in terms of C % rather than (7® , and 
denote the corresponding interference terms by G rather than . All the necessary Gfl’ 
and G-^ baie reac@ 

MO) _ r(2 - e) e 7 ® 

T?7 T(2 - 2e) ’ 

Gi? = jG™ = -ir(l + e )e’*G£ ) , 

5 Iii particular, for the function given in Eq. (13) of Ref. [55] . we have lim mc _^o h. 2 ? (5 = 1) = ^ — | 7 r 2 . 

6 Gi 7 differ from Gi 7 only for i = 3,4, 5, 6 . 

8 





















<7? = 5 G™ = 4(5-36- e 2 ) G™, 


t(0) 


(0) 


G. 


(l)bare 

27 


= -6 G 


(l)bare 

17 


92 


81 e 243 


1978 777tt 2 

+ — 


27185 


729 


e + 0(e 2 ), 


G 


(l)bare 

47 


1 Q/^'7/1 

—1 + —— + 43.76456245573869 + 94.9884724116 e 
3 e 2 243 e 


G 


(l)bare 

77 


+ KUi 

+ n b 
4 


16 447r 2 

+ 


612 


243 


243 


+ ni 


16 


81e 243 


4 264vr 2 

+ - 


2186 


729 


16 

81e 


+ 0.04680853247986 + 0.3194493123 e + Q(e 2 ), 


G, 


(l)bare 

78 


124 
3e + IT 

16 280 
9e + 


16 2 , 

( 212 

58 9 

64 

— 7T 2 + 

- — 

— 7T- - 

-C.'i 

9 

l 3 

9 

3 

16 2 1 

/ 382 

16 2 

160 

— 7T 2 + 

27 1 

hr- 

— 7T 2 - 

9 

9< 


G e+(D(e 2 ), 


A(l)bare 

^7(12) 


c A(l)bare 
-U Cr 7 ( 11 ) 


2096 39832 2x 

+ -^e + 0(e 2 ). 


- , _ (2.4) 

81 243 v ’ 

The last line of the above equation describes contributions from the so-called evanescent oper¬ 
ators that vanish in four spacetime dimensions 


Qu = (s L r r f n'y^v 3 Ta c L )(c L 'y fll Y 2 l fl 3 T a bL) - I6Q1, 

Q 12 = (sL7w7/ t2 7/ i 3 c i)(cL7 w 7 /i2 7 M ^) - !6Q2- 


In 


Gj}^ 0 , the three-particle-cut contributions alone (6 —> s^g) read 


G 


( 1)3 P 


27 


= -6 G 


( 1)3 P 


17 


4 

27 


106 

~81 


e + 0 (e 2 ). 


(2.5) 


( 2 . 6 ) 


In addition, several interferences need to be calculated with the 6 -quark propagators squared, 
to account for the renormalization of m b . We find 

1 21 + 4tt 2 

3^“ 


G 


(l)m _ 


27 


= -6 G 


(l)m _ 
17 — 


81 e 


+ 


1085 


81 


161 2 40 

-7T-(3 

972 27 s 


/59071 _ 1645^ _ 65 T_ , 

486 2916 81 81 1 K J 


{ 


G 


(0 )m 
47 


4 _ 50 - 2vr 2 94 - 3vr 2 - 32Cs _ 2 


3l + 2 + 


e + 


e 2 + 0{e A ). 


(2.7) 


9 9 

Our conventions for their global normalization will become clear through the way they enter 
the renormalized NNLO expression in Eq. (12.10)1 below. 

Some of the diagrams with Q\ insertions contain 6 -quark tadpoles that are the only source 
of 1/e 2 terms in &^ haxe ^ and 1/e terms in G^ m . Such divergences are actually necessary to 
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renormalize the 1/e 3 poles in Eq. (12.31) . These tadpole diagrams have been skipped in the NLO 
calculation of Ref. [33] because they give no contribution to the renormalized G ^ 7 . i.e. they 
cancel out after renormalization of rrif,. 

Among all the bare interferences given in this section, not only the NNLO ones are entirely 
new, but also G^l ^ e , G^ m and • The remaining LO and NLO results are extensions of 
the known ones by another power of e, as necessary for the current calculation^ 


2.2 Renormalization 

Our results in the previous subsection contain no loop corrections on external legs in the 
interfered amplitudes. Such corrections are taken into account below, with the help of on-shell 
renormalization constants for the 6 -quark, s-quark and gluon fields 

Z? s = 1 -jS,i> e ^‘r(£)^| + o(5 s 2 ), 


z° s = 1 + o(sa 


Zq S = 1 - n b «s s e e 7e T(e) + O(a^), 


( 2 . 8 ) 


where cu = t 2 - = 


and s = ^ 2 -e 7 . The QCD coupling g s and the Wilson coefficients Ci 

o 

are renormalized in the MS scheme: g bare = Z g g s , and C bare = JA CjZji. The corresponding 
MS renormalization constants can be taken over from the literature (see, e.g., Refs. HZIHSI) 


J 9 

Z\\ 

Z\2 

Z\z 

Z\A 

Z 15 

Z\§ 

Zyj 

ZlS 


z„ = 


1 + t (-T + I) + °(5. 2 ). 

1 - 2J t + 0(5 2 ), 


4 a. 


3e + ®( a s)i 

10 353 ) + 0(5 3 ), 


a s ( 


81 e * 1 243 e; 

-|Z 24 + oil (2^2 - , 

a l ( _ 8T? + 486^) + °( a s)» 


= a 


( 216^2 - eib) + °( a s)> 


216 e 2 

_ 1 7 1 ~2/ 22 332 \ 

6 Z 2 7 + Cl s ( gle 2 243e) ’ 

167 + 0 ( 5 S 2 ), 

+ ^( Q s)i 


648 e 
5 5 s 

12 e 


^ 1 ( 12 ) — ifr + &( a s)> 


z 77 

Z 21 

Z 22 

Z 23 

Z 24 

Z 25 

Z 26 

z 27 

Z 28 
z 2(11) 
■^ 2 ( 12 ) 


= 1 + 


16 a s 
3 e 


+ <?(5 2 ), 


6 -f + 0(a*), 

1 + 0(al), 


a 2 s (- 

2 5, 

3 e 


20 


27 e 2 


- + C^a 3 ), 


+ s s 2 + if) + C»( 5 3 ), 

a s (rh 1 + + ®( a s)> 

a s (36? + Toll) + °( a s)’ 

116 5s , ~2 ( -3556+744/ , 13610-44/ 1 , n( ~3\ 

81 e +U s l 2187 c 2 ~ l ” 2187 e 

^ + ^(S S 2 ), 

% + 0(5 s 2 ), 

(9(S s 2 ), 


(2.9) 


where f — rg + rib here, as we have skipped all the charm loops on the gluon lines. For the 
6 -quark mass renormalization, we use the on-shell scheme everywhere ( Z ° s = Z° s + O(al)), 
to get the overall mf pole in Eq. (11. 8 j) . 


t-, , • /v(0)bare /k(l)bare 

Exceptions are & 77 G* 


and Gto baie , for which sufficiently many terms in the e expansions have 


T 77 CU1U U 7 g 

been already found in Refs. [25ll27U37j . Our results agree with theirs, barring different conventions for the global 

1 + 0 (e) normalization factor (see the end of subsection 
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With all the necessary ingredients at hand, we can now write an explicit formula for the 
renormalized interference terms up to the NNLO (i = 1,°^ 


S,G“+S, 2 Gg> = Z° s Z° s Z 77 \5*s 3 ’Gy b ‘” + (Z° s ~l)s‘ 


>1 g'? 


Z/ z 4 (jr 47 a s S Lr i7 


+ a s (Z° s - 1) S 2e + Z„ Z° s 

,(0) , ~ e ry2 /A(l)bare 


a $ + a s s e G& )bare 


-I- rv 7 e e A(l)bare 
+ «s A8 S Ct 78 


E 

J=l,— j6j11j12 


Ejj s e 


G^ + a s s e Z 2 G ) 7 


+ ^( a s)l 


( 2 . 10 ) 


where G^J = 0 for j = 1,2,11,12. Once the above expression is expanded in a s , and 0( 3%) 
terms are neglected, all the l/e n poles cancel out as they should. Our final renormalized results 
at E 0 = m c = 0 read 


/"»(!) 

Or 27 — 

-6 G" 

= 

i 

A( 2 ) 

Ori 7 — 

1 A(2) 

6 G 27 

+ 

136 

27 

A( 2 ) 

^27 — 

/11792 

+ 

800 

V 729 

243 


243 


-In —, 

81 m b 


2 /i 94 + 87t 2 h 


In 2 — + 
m b 


9 

i 2 H— 
m b 


In — + 22.6049613485, 
rn b 


64 

729 


k Til 


2300 \ 

——m + 9.6604967166 n b ) In— - 14.0663747289 + 0.1644478609 nm 
243 J m b 


(54170 92 9 , 

W + 729" 1 - 1 ' 8324 ° 81161 "*' 


( 2 . 11 ) 


They are, of course, insensitive to conventions for the global 1 + O(e) normalization factor in 
Eqs. (I2.3p - fl2.7p . so long as it is the same in all these equations. In particular, it does not 
matter that our G 77 differs from the one in Ref. [25] by an overall factor of T( 1 + e) e 7£ . 

As already mentioned, the ni terms not marked by k in Eq. (12. lip agree with the previous 
calculations where both m c ^ 0 and m c = 0 were considered. In the case of the n b terms, 
the current result extends the published fit (Eq. (3.3) of Ref. [32] ) down to m c = 0. All the 
remaining terms are entirely new. 

8 Obviously, the renormalized G remain unchanged after replacing Z g —> Z gi Z t j —> Zjj and s —» 
on the r.h.s. of Eq. (12.101) and inside the on-shell constants (12.81) . 
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3 Impact of the NNLO corrections to (Q 7 , Q 12 ) interfer¬ 
ences on the branching ratio 

In the description of our phenomenological analysis, we shall strictly follow the notation of 
Ref. [34], where the relevant perturbative quantity 


p(e 0 ) = 

i,j =1 


(3.1) 


has been defined through 

T[b -» X^] Ei>Eo 
\V cb /V ub \ 2 r[6 -)• X*eV] 


v; s v tb 


v; 


c6 


6 a e 


7T 


W). 


(3.2) 


The relation between for i = 1, 2 and K i7 = a s + oi 2 Kff + 0(a g) is thus very simple 
qi) , ~2 t^( 2) , /o i ~3_ “sG-^ + asG-y +0(a s 3 ) 


a. 


+ 0(5. J ) = 


1 + 5 S (50 — 8tt 2 )/3 + £>(5 S 2 ) ’ 


(3.3) 


where the denominator comes from the NLO correction to the semilcptonic b —> X^ev decay 
rate. 

(2,} 

In the following, we shall write expressions for K \ 7 that are valid for arbitrary m c and E 0 
but incorporate information from our calculation in the previous section, where Eq = m c = 0 
has been assumed. For this purpose, four functions 

fNLo{z,S ) = Re r 2 1} (z) + 2 <j§j{z, 8 ), 
f q {z,5) = Rerf^) - § h$(z,5), 
f b (z ) ~ —1.836 + 2.608 z + 0.8271 z 2 — 2.441 2 In z, 

f c (z) ~ 9.099 + 13.20 2 ) — 19.68 z 2 + 25.71 z In z, (3.4) 

of z = m 2 /m 2 and 5=1 — 2 E 0 /m b are going to be useful. Explicit formulae for r£\z) and 
Rer^z) can be found in Eq. (3.1) of Ref. [43] and Eq. (26) of Ref. [30], respectively. For 
Jiqj {z, 6)1 we shall use a numerical fit from Eq. (13) of Ref. [39]. An analytical expression for 
027 \z,S) for Az < 1 — S (which is the phenomenologically relevant region) reads 

= -|-5(3 - 35 + 5 2 ) + \z5s s L s + 12 ~ ^ z 2 8 + 1 - 2 z)(s 0 L 0 - s s L s ) 

z ( 6 y o 

Ostt-2 _ 7 Q Q 

+ -^**(2 - 6) - ^z(6z 2 - 4z + 1 )(Lg - L 2 ) - |z5(2 - 5 - 4 z)L 2 , (3.5) 
with ss = \J (1 — 5)(1 — 5 — 4 z), s 0 = \/l — 4z, L b = In y/T ~ ? 0^~^ £ and L 0 = In 1 + ^0 4 ~ • 
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In the 5 = 1 case, <f>$ and h$ for z < \ are given by 


4*27 (a 1) — 
h?7(z,l) ^ 


2 12-87T 2 2 4 0 . r 2 tt 2 -7 8 2 , iNr2 4 2 3 

- 1 - z 2 H— z( 1 — 2z)s 0 L 0 H- -2 - z(6z 2 — 4 z + 1 )Li H—7r 2 z , 

27 9 3 ' ; 9 9 y 0 3 

— - -it 2 - 2.24 z 1/2 - 7.04 ^ + 23.72 z 3/2 + (-9.86 z + 31.28 z 2 ) In z. (3.6) 

27 9 v > \ j 


The functions f b (z) and / c (z) in Eq. (13. 4 p come from Eqs. (3.3) and (3.4) of Ref. [32] . 
respectively. These numerical fits (in the range z G [0.017,0.155]) describe contributions from 
three-loop b —* sy amplitudes with massive 5-quark and c-quark loops on the gluon lines. 

The ratio 0 = m 2 c /m 2 is defined in terms of the MS-renormalized charm quark mass at an 
arbitrary scale /i c . In practice, we shall use p c = 2.0 GeV as a central value. As far as the 
renormalization scheme for mb is concerned, we assume the following relation to the on-shell 
scheme 


^b,pole 

m b 


1 + oe s x m + 0(a 2 ). 


(3.7) 


In the IS and kinetic schemes, one finds x rn = |7rax and x m = 6 g(^ n ^1 + j , respectively. 
In our numerical analysis, the kinetic scheme is going to be used. 

Complete expressions for the NNLO quantities K jy and can now be written as follows 


K%\z,S) 


1 


S) + At + F l (z, 6) + I ^ - ~-K, 


( 94 




hi) 


6 


V 81 


U- 


34 


4‘ 78 / 27 bl 


= A 2 + F 2 (z, 5) - 3 f q (z, 5) + f b (z ) + fc(z ) + ^27 ( z > 7 lnz 


+ 


(8 L c 2x m ) z T (1 5)x m 

clz ad 


f N Lo(z,S ) + 




(i) 

27 


-A' 


(i) 


47 


208 

~81 


K& - —K. 
27 


35 


(i) 


254 

~81 


5948 

729 


J 2 

b'bi 


(3.8) 


where /3q 1 ~ 3 = 9, L b = In(p 2 /m 2 ) and L c = ln(^ 2 /m 2 ), while the relevant are collected 
in Appendix C. 

The expressions Ai + Fi(z,8) contain all the contributions that are not yet known for the 
measured value of m c . They correspond to those parts of the considered interference terms that 
are obtained by: (%) setting fi b = mb, /i c = m c and x m = 0, (n) removing the BLM-extended 
contributions from quark loops on the gluon lines and from b —>■ sqq'f decays (q = u,d,s ), 
except for those given in Fig. [2b. 

We define the constants A t by requiring that E)(0,1) = 0. Then we evaluate A t from 
Eq. (12. lip by setting there p = m b , n b = 0 and nrii = 3. Next, a replacement rq —> ui + ^/3q 1 = 7 
is done in the remaining rq-terms. Finally, Eq. (13.3p is used to hnd 


A x ~ 22.605, A 2 ~ 75.603. (3.9) 

These two numbers are the only outcome of our calculation in Section [2] that is going to be 
used in the phenomenological analysis below. 
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Apart from the condition Fj(0,1) = 0, everything that is known at the moment about the 
functions Fi(z,5 ) are their large -z asymptotic forms. They can be derived from the results of 
Ref. [35] @ Explicitly, we find 


70 

Fi(z,5) = — In 2 z + 


119 2 2 3 (i). . . , 

'W~9 n+ r™ ^ ) lnz 


493 


5 2 232 5 (1) . . 

- K H- (3 + -07S (5) 

2916 54 27 S 8^ 78 V ; 


F 2 (z,5) = 


A. + O 
4736 


729 


In 2 z + 


165385 1186 2 

+ -=^r7r 


+ J ( ln 2 + !)- 

- 2 ^£{s)~M + o n 


2187 729 

956435 2662 


It! + l Yl + \^ {S) + 


19683 2187 


7T + 


20060 

243 


C 3 “ 


1624 
243 ( 




(S) 


(3.10) 


The constant Y\ and the necessary 0 ^ functions are given in Appendices B and C, respectively. 
Let AB S1 , denote the contribution from Fi j2 (z,8) to B S1 . Then the relative effect is given by 


A B 


c ‘ 0) (w) F >M) + (<T'(w) - s c r M) fiM) 


S 7 


B 


~ U(z,5) 


1 


1 /m( 0 ) / 


s 7 


87T 2 


a 


(O)eff / 


7 (/h>) 


(3.11) 


For fi b = 2.0 GeV, we have a s (/i b ) ~ 0.293, C[°\/i b ) ~ —0.902, Cf\ii b ) ~ 1.073, and 
C , 7° )eff (^&) - -0.385. 

We shall estimate the contribution to B S1 that comes from the unknown U(z,S) by consid¬ 
ering an interpolation model where U(z, 1) is given by the following linear combination 


Lintei-p(^, 1 ) = X\ + X 2 f q {z,l) + ^X 3 + X 4 Z-^j f NLO (z,l). 


(3.12) 


The numbers aq are fixed by the condition 1/(0,1) = 0 as well as by the large-z behaviour 
of U(z,l) that follows from Eq. (13.101) . This determines x* in a unique manner, namely 
Xi — (—0.0502, 0.0328, 0.0373, 0.0309);. In Fig. 0] the function U- mteT p(z, 1) is plotted with a 
solid line, while the dashed line shows U asymp (z, 1), i.e. asymptotic large-z behaviour of the true 
U(z, 1). Note that yfz = m c /m b rather than z is used on the horizontal axis. The vertical line 
corresponds to the measured value of this mass ratio. The plot involves some extra approxima¬ 
tion in the region between yfz — 0.4 and \fz ~ 0.8 where we need to interpolate between the 
known small -z and large-z expansions of Ren, (z) (see Fig. 1 of Ref. EU). 

In Refs. [341142] the uncertainty in B S1 due to unknown m c -dependence of the NNLO correc¬ 
tions has been estimated at the ±3% level. The size of the interpolated contribution in Fig. [4] 

9 We supplement them now with the previously omitted large-m c contributions from the diagrams in Fig. 1 
in Ref. [35, or, equivalently, Fig. [3jo in the present paper. The effect of such a modification is numerically very 
small. 
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Figure 4: The interpolating function defined in Eq. (13.121) (solid line) and asymptotic behaviour 
of the true function U(z, 1) for m c 7Ub/2 (dashed line). The vertical line corresponds to the 
measured value of m c /mb- 


implies that no reduction of this uncertainty is possible at the moment. One might wonder 
whether the uncertainty should not be enlarged. Our choice here is to leave it unchanged, for 
the following reasons: 

(i) Our choice of functions for the linear combination in Eq. (13.12ft is dictated by the fact, that 
these very functions determine the dependence on z of the known parts of K\ 7 and K 27 . 
The known parts are either those related to renormalization of the Wilson coefficients 
and quark masses (in the terms proportional to Lb and L c ) or the renormalization of 
a s (the function f q parametrizes the considered correction in the BLM approximation). 
It often happens in perturbation theory that higher-order corrections are dominated by 
renormalization effects. If this is the case here, the true U(z,l) should have a similar 
shape to U mterp (z, 1). 

(ii) The growth of U inteip (z, 1) for m c > m b /2 is perfectly understandable. In this region, 
logarithms of 0 from Eq. (13.10p combine with Lb from Eq. (I3.8[) . and the asymptotic large- 
m c behaviour of K (] 2 ^ 7 is determined by In (n b /m c ) and ln(// c /m c ) only (see Eqs. (5.12) and 
(5.14) of Ref. [35]). Thus, the growth of the correction for large z can be compensated by 
an appropriate choice of the renormalization scales, which means (not surprisingly) that 
the dangerous large logarithms can get resummed using renormalization group evolution 
of the Wilson coefficients, masses and a s . 

(iii) Our ±3% uncertainty is going to be combined in quadrature with the other ones, which 
means that it should be treated as a “theoretical lcr error”. To gain higher confidence 
levels, it would need to be enlarged. 

(iv) In the considered interference terms id 17 and K 2 7 , the dependence on 5 is very weak in the 
whole range 6 G [0,1], both at the NLO and in the BLM approximation for the NNLO 
corrections. Specifically, changing S from 1 (Eq = 0) to 0.295 (Eq = 1.6 GeV) results in 
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1 2 

3 

4 

5 

6 

7 

8 

9 

10 

total 

-0.6% +1.0% 

-0.2% 

+2.0% 

+1.0% 

+1.6% 

+2.1% 

-0.5% 

+0.2% 

-0.4% 

+6.4% 


Table 2: Shifts in the central value of B S1 for E 0 = 1.6 GeV at each step (see the text). 

modifications of }nlo by + 0 . 2 % and f q by + 1 . 0 %, respectively, for the measured value 
of m c . The corresponding changes at m c = 0 amount to —0.7% and —2.4% only. Thus, 
our estimates made for 8 = 1 are likely to be valid for arbitrary 5. 

( 2 ) ( 2 ) 

In the phenomenological analysis below, we shall take K\ 7 and K, 27 as they stand in 
Eq. (BSD , replace the unknown Fi(z,5) by Fj inteip (^, 1) interpolated analogously to Eq. (13.12j) 

= -23.75 + §/,(,, 1) + gW*.!), 

Ef terp (+1) = -3.01 - ^f q {z,l) + (-10.34 - 9.55 z-^j f NLO (z, 1), (3.13) 

and include a ±3% uncertainty in the branching ratio due to such an approximation. 

4 Evaluation of B S1 in the SM 

In the present section, we include all the other corrections to B S1 that have been evaluated 
after the analysis in Refs. [341142] , Next, we update the SM prediction. To provide information 
on sizes of the subsequent corrections, the description is split into steps, and the corresponding 
modifications in the branching ratio central value are summarized in Tab. [2] The steps are as 
follows: 

1. We begin with performing the calculation precisely as it was described in Ref. [34] but 
only shifting from B(B -+ A" s y) to B SJ , which amounts to CP-averaging the perturbative 
decay widths. No directly CP-violating non-perturbative corrections to B(B -+ Xg'y) were 
considered in Ref. [34]. It was not equivalent to neglecting them but rather to assuming 
that they have vanishing central values. A dedicated analysis in Ref. [48] leads to an 
estimate of 0.4 ± 1.7% for such effects. 

2. The input parameters are updated as outlined in Appendix D. In particular, we use results 
of the very recent kinetic-scheme fit to the semileptonic B decay data [49j . 

3. Central values of the renormalization scales (/i c ,/q,) are shifted from (1.5, 2.5) GeV to 
(2, 2) GeV. Both scales are then varied in the ranges [1.25, 5] GeV to estimate the higher- 
order uncertainty. In the resulting range of B S1 , the value corresponding to the (2, 2) GeV 
renormalization scales is more centrally located than the (1.5, 2.5) GeV one, after per¬ 
forming all the updates 1-10 here. It is the main reason for shifting the default scales. 
The (2, 2) GeV choice is also simpler (both scales are equal), and /i c is exactly as in the 
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Figure 5: Interpolation of p^ 2)rem in m c as in Fig. 2 of Ref. [34J but with updated input parameters 
and with renormalization scales shifted to (// c , fi b ) = (2, 2) GeV. In addition, the thick solid (red) 
line shows the case with the presently known boundary condition at m c = 0 imposed. 


fit from which we take m c (/t c ) (Appendix D). As far as /if, is concerned, it should be of 
the same order as the energy transferred to the partonic system after the 6 -quark decay. 
For the leading b — > sy contribution from the photonic dipole operator P 7 , this energy 
equals to \m b which gives 2.3 GeV when one substitutes m b = m bt kin from Appendix D^j 
Rounding 2.3 to either 2.5 or 2.0 for the default value is equally fine, given that the 
observed /^-dependence of B S1 is weak (see Fig. [ 6 ]), and our range for /if, is [1.25, 5] GeV. 

4. In the interpolation of p^ lcm (see Ref. |34j for its definition), we shift to the so-called 
case (c) where the interpolated quantity at m c = 0 was given by the (Q 7 , Q 7 ) interference 
alone. 

5. The m c = 0 boundary for p^ 2>rc,u j s updated to include all the relevant interferences, 
especially the ones evaluated in Section [2j The thick solid (red) line in Fig. [5] shows the 
new p^ rem in such a case, while the remaining lines are as in Fig. 2 of Ref. [31] (somewhat 
shifted due to the parameter and scale modifications only). 

6 . At this point, we abandon the approach with m c -interpolation applied to the whole non- 

BLM correction p ^ rem . As before, the penguin operators Q^,...,6 and the CKM-suppressed 

( 2 ) ( 2 ) 

ones Q 12 are neglected at the NNLO level. The corrections K\ 7 and K 27 are treated 

( 2 ) 

as summarized at the end of the previous section. For K) 8 , the complete results from 
Refs. |36|[37] are included. K 77 is made complete by taking into account its exact m c - 
dependence [29] 150]. in addition to the previously included terms. For the NNLO inter¬ 
ferences among Q 1 , Q 2 and Qg, only the two-body final state contributions are present at 
this step. They arc infrared-finite by themselves, and given by products of the well-known 
NLO amplitudes r [ 1 (see Eq. (3.1) of Ref. [T3]) whose imaginary parts matter here, too. 

10 The measured photon spectra are also peaked at around 2.3 GeV, which confirms the leading role of the 
two-body partonic mode. 
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7. Three- and four-body final state contributions to the NNLO interferences among Q i, Q 2 
and Q 8 are included in the BLM approximation, using the results of Refs. [3TII38H39] . Non- 
BLM corrections to these interferences remain neglected. The corresponding uncertainty 
is going to be absorbed below into the overall ±3% perturbative one. 

8 . Four-loop <5i,...,6 —>> Qs anomalous dimensions from Ref. [19] are included in the renor¬ 
malization group equations. 

9. The LO and NLO contributions from four body final states are included [231124] . They are 

not yet formally complete, but the only neglected terms are the NLO ones that undergo 
double (quadratic) suppression either by the small Wilson coefficients ^ or by the 
small CKM element ratio \ V* s V u t,\ / The uncertainty that results from neglecting 

such terms is below a permillc in B S1 . As far as the CKM-suppressed two-body and 
three-body contributions are concerned, the two-body NLO one has already been taken 
into account in Ref. [BUJ. The remaining NLO and NNLO ones (also those with double 
CKM suppression) are included at the present step. Their contribution to B S1 is below 
a permille. However, the branching ratio B^ [51] receives around 2% enhancement from 
them. 

10. We update our treatment of non-perturbative corrections. The O (a s A 2 /m 2 ) correction to 
the ( Qy , Qt) interference from Ref. [40] replaces the previous approximate expression from 
Ref. |52j. Moreover, we include a similar correction [4T|]53] to the charmless scmilcptonic 
rate that is used for normalization in [P(i%) + N(Eo)] (see Eqs. (ID. 2D and (ID.411 in 
Appendix D). In consequence, the previous (tiny) effect in N(E 0 ) gets reduced by a factor 
of around 4. Finally, our treatment of non-perturbative effects in interferences other than 
gets modified according to Ref. [14] . A vanishing contribution to the branching 
ratio central value from such corrections is assumed, except for the leading 0(y\ 2 /m 2 c ) 
one 0 where m c is fixed to 1.131 GeV. At the same time, a ±5% non-perturbative 
uncertainty in the branching ratio is assumed, as obtained in Sec. 7.4 of Ref. CHI by 
adding the relevant three uncertainties in a linear manner FH 

Our final result reads 

B™ = (3.36 ± 0.23) x 10 “ 4 (4.1) 

for E 0 = 1.6 GeV, where four types of uncertainties have been combined in quadrature: 
±5% non-perturbative (step 10 above), ±3% from our interpolation of F lt 2 (z,S) (Section |3]), 
±2.0% parametric (Appendix D), as well as ±3% from higher-order perturbative effects. 

The latter uncertainty is assumed to account for approximations made at the NLO and 
NNLO levels, too. In the NLO case, it refers to the doubly suppressed terms mentioned in 
step 9 above. I 11 the NNLO case, it refers to neglecting the penguin operators at this level, 
and using the BLM approximation in step 7 above. If we relied just on the renormalization- 
scale dependence in Fig. [6] (with 1.25 GeV < /i c , ^ < 5 GeV), we could reduce this uncertainty 

11 If their ranges were treated as ler ones and combined in quadrature, the uncertainty would go down to 
3.3%. 
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Figure 6: Renormalization scale dependence of B S1 in units 10 -4 at the LO (dotted lines), NLO 
(dashed lines) and NNLO (solid lines). The upper-left, upper-right and lower plots describe the 
dependence on /r c , /p, and /x 0 [GeV], respectively. When one of the scales is varied, the remaining 
ones are set to their default values. 





to around ±2.4%. However, apart from the scale-dependence, one needs to study how the 
perturbation series behaves, which is hard to judge before learning the actual contributions from 
Fi >2 (z, 5). Thus, we leave the higher-order uncertainty unchanged with respect to Refs. (341142] , 
Our treatment of the electroweak corrections [55] remains unchanged, too. 

The central value in Eq. (14.11) is about 6.4% higher than the previous estimate of 3.15 x 10 -4 
in Refs. [54M2] . Around half of this effect comes from improving the m c -interpolation. As seen 
in Fig. O the currently known m c = 0 boundary for the thick line is close to the edge of the 
previously assumed range between the curves (a) and (b). It is consistent with the fact that the 
corrections in steps 4 and 5 sum up to 3% being the previous “la” interpolation uncertainty. 
The m c = 0 boundary has been the main worry in the past because estimating the range for 
its location was based on quite arbitrary assumptions. It is precisely the reason why no update 
of the SM prediction seemed to make sense until now, given moderate sizes of the other new 
corrections. 
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5 Conclusions 


We evaluated 0(a 2 s ) contributions to the perturbative T(b —* X s ^) decay rate that originate 
from the (Q 7 , Qip) interference for m c — E (} — 0. The calculation involved 163 four-loop massive 
on-shell propagator master integrals with unitarity cuts. Our updated prediction for the CP- 
and isospin-averaged branching ratio in the SM reads = (3.36 ± 0.23) x 10 -4 . It includes 
all the perturbative and non-perturbative contributions that have been calculated to date. It 
agrees very well with the current experimental world average £>“ p = (3.43 ±0.21 ±0.07) x 10~ 4 . 
An extension of our analysis to the case of B^ and an update of bounds on the Two Higgs 
Doublet Model is going to be presented in a parallel article [5T] . 
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Appendix A: Massless master integrals 

In the course of this work, it has been necessary to compute a number of massless scalar integrals 
with various unitarity cuts. All of them are depicted in Fig. [7] and Tab. [3] They occur after 
applying the large mass expansion for pi <C m 2 , as well as in the decay rate calculation itself. 
Apart from the four-loop diagrams with four-particle cuts, and the four-loop diagrams 4L3C1, 
4L3C2 and 4L3C3 with three-particle cuts, values of all our master integrals can either be 
found in the literature [56H60] or obtained using standard techniques described, for instance, in 
Ref. |64j. Let us note that the results for all the massless propagator four-loop master integrals 
in Refs. [EH 66] are not sufficient here because they correspond to sums over all the possible 
cuts, while certain cuts need to be discarded in our case. 

I 11 the following, we explain our computation of the four-particlc-cut master integrals in 
dimensional regularization with D = 4 — 2e. The total momentum is q = pi +P 2 + P 3 +Pa, and 
we have p 2 = 0 for i = 1,... ,4. Moreover, all the internal lines are massless. The momenta 
are in Minkowski space, and we tacitly assume that all the propagators below contain an 
infinitesimal +ir] with rj > 0. We also define the invariants 

Sijk... = (Pi + Pj+Pk + • • -) 2 • A- 1 ) 
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Figure 7: The massless four-particle-cut diagrams calculated in the course of this work. 


We therefore have s 42 + s 43 + 6’ 14 + s 23 + s 24 + s 34 = q 2 as a constraint from overall momentum 
conservation. 

Our convention for the loop measure is 



and we define the prefactor 

c 1 

r (An) D/ 2 T(l-e) ' 


(A.2) 


(A.3) 


Note that our definition of Sr is different from the one in Eq. (4.13) of Ref. [57] . 

As far as integration over the four-particle massless phase space in D — 4 — 2e dimensions 
is concerned, we closely follow Ref. The phase space measure reads 


dPS A = 


d 0 - 1 


Pi 


(2U d -' 2 Ek ■ ■ ■ (2ir)D 5lD,( " “ Pl “ P2 “ Ps - Pi) ■ 

It can be rewritten in terms of invariants and angular variables according to 


\4-3D ( 1— if ol-fr 


2 4 2 (—A 4 ) 2 0(—A 4 ) dVtr)_2 dVtn_^ 


(A.4) 


dPS A = (2vr) 4 - 3 " (q 2 y 

xd(q~ — S12 — s 43 — s 44 — s 23 — s 24 — s 34 ) ds\2 dsi 3 ds\± ds 23 ds 24 ds 34 , 


(A.5) 


with the Gram determinant 


A 4 = A(si 2 s 34 , si 3 s 24 , si 4 s 23 ), X(x, y, z) = x 2 + y 2 + z 2 - 2 xy - 2 xz - 2 yz . (A.6) 


21 























Table 3: The massless two- and three-particle-cut diagrams used in the course of this work. 

It turns out that integration over angular variables is trivial in all the cases we encounter 
here, and we can use 

/' 2v t d / 2 

r B= wm- (A7) 

Performing the angular integration, and furthermore applying the steps explained in Ref. [57j 
to factorize the phase space measure, we arrive at 

dPS 4 = - —- 35 --—-—-- dt dv d\ dzi dy lu dy vlM 5(1 - y X2 u) (A. 8 ) 

(4tt) —(1 - 2e)T(l - e)T 2 (l - e) 

t~ e (i - 1)- £ u- (i - u)' e x'^ e (i - x)“*“ e *r (i - y- 2 t ylu e (i - i/m) 1 " 26 . 
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All the integration variables t, v, x, z i, y 134, and 1/1234 run from 0 ... 1 and originate from 


_ 2 

$ijk... Q Vijk... 1 

1/13 — (l/l3,6 — 1/13, a)X + 1/13, a , 

yi 2 = y 134 zi t , 

1/13, 6/a = B ±y/B 2 - C , 

1/23 — 1/134 O , 

B = Visa {tv + vtz 1 ) , 

1/14 = 1/134 0 ^ , 

<7 = l/?34(tU - Vt2i) 2 , 

1/24 = 1/134 £l t , 

y/B 2 - C = 2 1/134 VtVt Vv y/v y/z[ 

1/124 = 2l (1 - Vl3AV) , 

1/13,6 - 1/13, a = 2\AB 2 - O , 


where t — 1 — t, and analogously for all the other variables. The substitutions (1A.9I) should be 
done in the integrands, too. 


A.l Results for the four-particle-cut master integrals 

We are now in position to present results for the four-particle-cut diagrams depicted in Fig. [71 
Normalization factors are extracted according to 

W* = 2-7T e i7re (g 2 ) ai " 4e 7 4L4 « , (A.10) 

where the a* follow from dimensional considerations. One finds a* = (2, 2,1, —1, 0, —1, —1, 0)* 
for i — 1 ,..., 8. 


We start with I^laci-, 


Ialaci 


which yields 
Ialaci 






, r(e)r 2 (i 
(47 t) d/2 r(2 - 2e) 


(trfdi 


r(e)r 9 (l - e)r(l - 2e)r(2 - 3e) 
T 2 (2 — 2e)r(3 — 4e)T(4 — 5e) 


The next integral to consider is Ialac 2 , 


Ialac2 


and we get 
IALAC2 






e ine T(e)T 2 (l 

(4tt) d/2 T(2 - 2e) 


WT 


T(e)T 10 (l — e)T(2 — 3e) 

T 2 (2 — 2e)T(3 — 3e)T(4 — 5e) ' 


(A.ll) 


(A-12) 


(A.13) 


(A.14) 
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We proceed with Ialac 3 , 


ItLaC3 JdPSaJldk} (7,; - /(| + /J; , - p 4 )2 ( Pl + p 2 + p 4 p 

= (?r(2-~2 (A.15) 

and arrive at 

Wa = F rip °-2€)r(4 ( -5t ) €) 3 ^( 1 ,1 - 2 - 3t; 2 - 2t, 4 - 5e; 1). (A.16) 

The expansion of Ialac 3 in e is conveniently done with the package HypExp [5B:,iS9], 


UL4C3 — 


4e 


+ 


37 


vr 2 \ 
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/198241 
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192175tt 2 
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12995?r 2 3521vr 4 


35tt 2 „ \ (13677 
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253vr 2 29 t r 4 
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175197T 4 

960 
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1481tt 6 
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V 32 


1287 67 9> 
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16 

315 


Cs 


e 3 + 


- 71Cs 

144 S 

/2597477 


19139 


925 


C 3 + -r^Cs + 17 °C 
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V 128 
2049C 5 


(A.17) 


We now move to Lilac a, 
Ialaca = / dP S 4 


[dk\ 


k 2 (k + pi + p 3 + p A ) 2 (pi + p 3 ) 2 (pi + p 2 + Pa ) 2 (pi + P2) 2 


e^ £ r(e)r 2 (l -e) 
(4i) D/2 r(2 - 2e) 


(A 


2 \ -3-« 


dPSA Ui 3 a Ul24 Vl3 Ul2 j 


(A.18) 


which does not reveal a closed form since we cannot avoid y 13 in the integrand. We therefore 
compute it from the following two-fold Mellin-Barnes representation 


C^+lOO C9+4 00 


; 2+ic 


/4L4C4 


r(e)T 6 (l — e)T(—e)T(l — 3e) f dz 3 f dz 2 


2 ni / 27 Vi + Z2 ~ e )^( _e ~ Zl ~ 2 2 )r(zi) 


x 


T(—2e)T 2 (2 — 2e) 

c-^—i 00 C2 —i 00 

r(i — *i)r(i — 2e — ^1) r(-* 2 )r(i + ~2)r(— 1 — & — Z2)r(i — e + z 2 ) 


T(2 — 5e — zi) 


T(1 - 3e + z 2 ) r(— e - z 2 ) 


(A.19) 


The integration contours in the complex plane can be chosen as straight lines parallel to the 
imaginary axis. The integral is then regulated [67] for ci = 1/2, C 2 = —1/4, and e = —7/4. We 
perform an analytic continuation to e = 0 with the package MB.in [67], which is also used for 
numerical cross checks. The expansion of Ialaca in e reads 


Ialaca 


1 1 
W + V + 




13tt 2 33 \ 1 ( 13tt 2 397tt 4 

—--2 c #H 20 -—-W 
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—66^3 I —hi 48 —-—-— 1980 H-vr C3 —-C5 

s 7 e V 3 360 4 s3 2 s5 

130vr 2 397 t r 4 24539 tt 6 „ 2a 897 a2 _ . 

112 “ — - W - 5286 + 13l2r 6 + T$ - 1374& I £ 


0(e 2 


(A.20) 


The next integral, In 4 Cb , with 

W 5 = [dPsJ[dk] 


k 2 (k + p 4 ) 2 (k+ p 4 +p 2 + p 4 ) 2 (p 2 +Pi) 2 


(A.21) 


(2) - 2 - e r r ix 

(4vr) D/2 r(l-2e) V ' J J 


[2/12 + 07 2/14 + o: 2/24] 1+e 2/23 ’ 


can again be expressed to all orders in e. One first integrates over x, and finally finds 


* 4L4C5 


r(e)r 6 (i-e)r 3 (- e ) r r(i-e) 
T(2 — 5e)T(2 — 2e) [r(2 - 2e) 
T(1 — 3e) 


3 F 2 (l,l-e,l-*2e; l + €,2-2e; 1) 


(1 - 3e)r(l - 4e) 

The expansion of hL 4 C 5 in e reads 


F 2 ( 1,1 — e, 1 — 3e ; 1 + e, 2 — 3e ; 1) 


. (A.22) 




*(««•* «, 

Also the next integral, / 4 l 4 c 6> with 

/4L4C6 = [dPsJldk] 


403vr 4 140 2 


7997T 6 

^C3 + 798 C5 + w 


125< 3 2 ) £ + Off 2 ) . (A.23) 


k 2 (k - p 2 ) 2 (k + p 4 ) 2 (k + pi + p 4 ) 2 (pi + p 2 ) 2 
e i7re T(2 + e)T 2 (—e) / 2 \-3-e ~ 1 


(A.24) 


(4ir) D/2 r(-2£) 


■(«*) 


cTPSit dx dy 


0 0 


® I/24 + 2/ 2/14 + xy yu} 2+e 2/12 ’ 


reveals a closed form which, however, turns out to be more complicated. One first integrates 
over x and y , and finally finds 


'4L4C6 


T(e)T 6 (l - e)T 2 (—e)T(—1 - 3e) 


- T(1 - 2e)T(e) - 2 T 2 ( 1 - 2e)T(2e)T(l + e) 


T(1 — 5e)T(2 — 2e)T(l —4e) 

-2T(1 - 2e)T(l + e) (^ (0) (1 - e) - ^ (0) (e) - ^ (0) (1 - 4e) + 2^ (0) (1 - 2e) + 7) 

4T 2 (—2e) 


—4T(—e) 3^2(1, — ■e, — e ; 1 + e, 1 — e ; 1) — 


T(—3e) 


yF 2 {-e, -e, -e; -3e, 1 - e ; 1) 
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r 2 (l-e)r(l-4e) 

(l + e) 2 T(l -3e)r(-2e) 
r 2 (l -2e)r(l + e) 
F(W) 4Fs( 


1 ^ 3 ( 1 ,1 — e, 1 — e, 1 + e ; 2 + e, 2 + e, 1 — 3e ; 1) 


4 -^ 3 (1) 1,1 — 2e, 1 — 2e ; 2,2,1 — 4e ; 1) , 


(A.25) 


where ip(°\z) = ^lnT(z). The expansion of I 4 lac& hi e reads 

~ 5 5 /35 79tt 2 \ 1 ( 65 79tt 2 \ 1 /275 553vr 2 

We - + ^ 36 

6437T 4 \ 1 / 665 1027vr 2 643vr 4 , 1301 „ 2590 \ 

+■mo + 58C d 7 + {-— + — -mo- 406C3 + —" 6 - — Cs ) 

f 2315 43457T 2 4501vr 4 63229 tt 6 „ _ 1301 2a , >2 2590 \ 

+ ( —-— + -Jieo - + 272I6O + 754(3 “ ~ 1 * + 1884( 3 + ~ Cs ) £ 


(A.26) 


'—(<2 'i 

The next integral, Ialacti has not been necessary for the actual calculation of G\ 7 and G 27 
because it stems from diagrams where the charm quark loop is cut. However, we still give the 
result, as it is the most complicated integral, and might be useful for future computations of 
other interferences. The difficulty is due to the fact that one cannot avoid 2/13 in the integrand, 
and the resulting Mellin-Barnes representation is four-dimensional. Starting from 


= dPS 4 / [dk] 


k 2 (.k - pi) 2 (k + p 2 +p 3 + Pa) 2 {k + p 3 + Pa) 2 (pi + P2 + P3) 2 


e l7Te T(2 + e)T^(—e) , 2 \- 3 -e 
(4ir) D/ 2 r(-26) (q) 


x / dPS 4 dx dy 


[ 2/34 + X ( 2/13 + yu) + y (P 23 + 2 / 24 ) + xy 2 / 12 ] 2+e 2/123 ’ 


(A.27) 


we first integrate over x and y, and find the following Mellin-Barnes representation. 


Ulaci — 


00 C2+ioo 

T(e)r 5 (l — e)T(—e)T(l — 3e) f dz x j dz 2 

T(—2e)T(2 — 2e) J 2n~i J 


q-ioo C2 - 


T(—e — zi)T(l — e + 2 i)T( 1 — 3e + z 4 — £ 2 )T(1 — 2e — z 2 )r(—e — z\ + z 2 ) 
r(l — z 2 — 3e)T(l — z 2 — 4e)T(l + z 2 — e) 

T(- Zl )T(l + zi)T(-g 2 )T(l + z 2 )T(-z 2 - e)T(^ 2 - e) 

T(1 — z\ — 3e)T(2 + z\ — 3e) 

ci+zoo C2+ioo 03+100 

2T(e)T 5 (l — e)T(—e)T 2 (l — 3e) f dz x f dz 2 j dz 3 

' T(1 - 5e)T(l - 2e)T(—2e)T(2 - 2e) J J 2n~i J 2rd 

q-ioo C2 —i 00 03—100 
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x r(—^i)r(i + z\ — z 3 )r(-2 2 )r(i + z 2 )v(—z\ + z 3 — P)P{z 2 — e)r(—^2 — ^3 — e ) 

r(l — z\ + z 2 + z 3 — 4e)T(2 + — z 3 — 3e) 

r(^ 3 )r(i — 4e + z 2 + z 3 )r(i — 2e + 2i)r(—^ + ^2 + ^3 — e )r(i — e + zi — *3) 

r(l + z 3 — 3e)r(l + z 2 — e) 

c^+ioo C2+ioo C3+4 oo 04+400 

r(e)r 5 (l — e)T(—e)r(l — 3e) f dz\ f dz2 f dz 3 f dz 4 
+ r(l - 5e)T(l - 2e)T(—2e)T(2 - 2e) ] ] 2^1 ] 

q-ioo C2 —i 00 eg—400 04 —400 

x r(-z 3 )r(z 3 — ^i)r(-z 2 )r(i + z 2 )v(— * 4 )r(i + z 4 + z 4 )r(z 2 — e) 

r(l - z 4 + z 2 + z 3 - z 4 - 4e)r(l + Zi — z 3 - e) 
r(l - 6 + zQTQsi - z 2 - z 3 - e)r(—^1 + z 3 - z 4 - e)T(-zi + Z2 + z 3 - z 4 -e) 

r(l + £3 — e)r(l + z 2 — e) 

xT( — Z\ — 2 e)T(l + 22 + ^3 — 2 e)r(l + z 4 — z 3 z 4 — e) . (A.28) 


The expansion of I 4 l 4 c 7 in e reads 
~ 2vr 4 


J4L4C7 


45e 


+ 


+ 2tt 2 C 3 - 58Cs^ 


+ 16^3 - 464C5 + 84CI - e + (9(e 2 


(A.29) 


We have also derived an alternative, seven-fold, Mellin-Barnes representation for I 4 l 4 c 7 and 
used it to confirm (IA.29|) numerically with the help of the code MB. m E2J. 


The last integral, I 4 l 4 cs , reads 
/4L4C8 = [dPS 4 f[dk] 


(4vr) D/2 r(l -2e) 


k 2 (k + p 1 + p 2 + p 4 ) 2 (k + pi+ p 2 ) 2 (pi +p 3 + p 4 ) 2 

1 

dPS 4 [ dx 


(A.30) 


e^r(l + e)T(-e)r(l-e) / 2 ,_ 2 -e 


V)' 


[yu + xy 14 + x y 24 ] 1+e 7/134 


Again, one first integrates over x, and finally finds an expression involving a one-dimensional 
Feynman parameter integral 


'4L4C8 
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T(1 - 3e)T 5 (l - e)T 4 (—e)T(l + e) 


2 T(2 — 5e)T(2 — 4e)T(2 — 2e) 
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2T(3 — 5e)T 2 (2 — 2e)T(—2e) 


3 + 2 ( 1,1 — e, 2e ; 1 + e, 1 + 2e ; 1) 
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x Jdt t 1 ~ 2e (1 -t)~ e 2 Fi(l,2-4e; 3 - 5e; t) 2 F 1 (l,l-e; 2-2 e;t) . (A.31) 
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The expansion of / 4L4C8 in e reads 
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A.2 Results for the three-particle-cut master integrals 

In this section, we describe onr computation of the three-particle-cut diagrams 4L3C1, 4L3C2 
and 4L3C3. Similarly to Eq. (lA.IOIh we extract the normalization factors according to 

hLza = 2vr e 2nie A 4 ( q 2 ) bi ~ 4e I 4L3Ci , (A.33) 

where the b\ again follow from dimensional considerations. One finds b\ — 0 and b 2 = —1. For 
4L3C3, we have used a different method, as explained below. 

The kinematics and the phase space measure are much simpler in the three-particle case, 
compared to the four-particle one. The total momentum is q = p 4 +p 2 + p 3 , and we have p 2 = 0 
for i — 1,..., 3. We define the invariants 


Sijk... = (Pi + Pj + Pk + • • -) 2 


(A.34) 


as before, and have Si 2 + s 33 + s 23 = q 2 as a constraint from overall momentum conservation. 
The phase space measure 
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(2n) D ~ 1 2E 1 (2tt) d " 1 2 E 3 
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(A.35) 


is again taken over from Ref. 0- After integration over angular variables one finds 


dPS 3 


27rS 2 T 2 (l -e) {q 2 f~ 2e 
T(2 - 2e) 


dy \2 dy l3 dy 23 y 12 e y 13 y 23 5(1 - y 12 - y V3 - y 23 ). 


The integration variables y 42 , y 43 , and y 23 run from 0... 1, and originate from Sij = q 2 yij. The 
latter substitutions have to be made in the integrands, as well. 


Our first three-particle-cut integral I 4 l 3 Ci reads 


Ialzci = / dPS 3 / [dk 4 \ / [dk 2 \ 
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It can be expressed in a closed form valid to all orders in e. One first integrates over x, and 
finally finds 
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T(2 — 4e)T(l — 3e)T(2 — 2e) 
r(l-2 e)r 5 (-e)r(2e)r 5 (l-e) 
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The expansion of I^lzci in e reads 
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The next three-particle-cut integral is Ial 3 C 2 , 
IaL 3C2 = IdPS, [ [dh] [ [ dk 2 ] 


{h +Pi+ P2) 2 k\ (hi -k 2 + Pi ) 2 {h - k 2 ) 2 (k 2 + p 2 ) 2 k\ ' 

(A.39) 


Despite the fact that p 3 does not appear in the integrand, the result of the integral is quite 
lengthy. In the end, we find the following expression that involves a one-dimensional Feynman 
parameter integral: 
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The expansion of I 4 L 3 C 2 in e reads 
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For the last integral I 4 L 3 C 3 , we employ a different approach. Due to the structure of the 
integrand, it is not possible to find a regulated Mellin-Barnes representation. Therefore, we 
begin with evaluating an integral Ial 3 C 3 ' defined as 

I4L3C3' — [dPS 3 [ [dki] I [dk 2 ] 
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0 0 

Again, we extract the normalization factor according to 
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The above quantity can be expressed in terms of a one-dimensional Feynman parameter integral 
as follows: 
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The expansion of Ial 3 C 3 ' in e reads 
~ 1 
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v 90720 

The original integral Ial 3 C 3 can then be obtained by relating it to 1^303' with the help of 
integration-by-parts identities. 


Appendix B: Relation to Ref. [43] 

Several decimal numbers in subsection 12.11 can be related to the quantities encountered in 
Ref. [12] as follows. In the finite part of G$ hare in Eq. (12.41) . we have 


43.76456245573869 = Y 1 = 


19039 11 

- -(- -7T 

486 27 


71 


16 1 

_ — X b -|—Re[a(l) — 26(1)], 

9x/3 27 6 6 1 U WJ ’ 


0.04680853247986 = Y 2 = 2Re6(l) 


243’ 


(B.l) 


where 


v 9 7T 2 2 1 l(1) fl 
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8 5 3 s3 RT V6 


„ 16 164 „ 16 3007T + 647T 3 32^^3-72 

Rea(l) =- 1 - 7 r-C3-7=-h 

v ; 3 405 9 ^ 


„ 320 632 

Re 5(1) =-1- 

v ' 81 1215 


71 


47r 
3^3 


135^3 

iWA, 

45 ^ Vey 


405 




and 


= ypum)- 


tp 


(B.2) 


(B.3) 


The above exact expressions for Xb and Rea(l) are new. They come from the three-fold 
Feynman parameter integrals in Eqs. (3.2) and (3.3) of Ref. [53]. 

In the 7-part of G^ bare in Eq. (12.31) . we have 

- 67.66077706444119 = -~Y X - 103762 + il^ 2 _ 

3 2187 27 27 


5.17409838118169 = 


2 11384 

- Y 2 + -. 

3 2187 


(B.4) 


Finally, in the coefficients multiplying ln(/x/m&) in Eq. (12. lift . we have 


1.0460332197 = 
9.6604967166 = 


■ 5 * 


37708 304 


-1-7T~ 

729 27 


4 7088 

— Y 2 T-. 

3 729 


(B.5) 
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Appendix C: 


NLO results of relevance for Section [3 


The NLO quantities that occur in Eq. (13. 8 ft are given by 

k£> = - 6 K$ = Rer? - ™L b + 2<f>$(8), 

k£> = RerP + ^L b + 2^(S), 

~ + f^ + 40ff(<5), 

= y-|;vr 2 + ^L 6 + 2 4V), 


(C.l) 


where and r[ l> can be found in Eq. (3.1) of Ref. [43]. The function 027 has been already 
given in Eq. (13,5j) here. The remaining ones read 




(5) 





2. 2 - 7. . 31 10. l- 2 2-, 1... ... - 

— hr 5-In 5-1-5 H— 5~ - 5 3 H— 8(8 — 4) In 8, 

3 3 9 3 3 9 3 v ' 

07 A W + 07 B W, (C.2) 

( 3%/3 ~ a ) + s' 53 ~ ts 152+ h s+ 1 o +26+3) arcta1 ' 2 vfrI 
l o + 43 + VjtI aict<m VitI 

345 2 + 595 - 18 5 2 In 8 , 433 5 3 + 429 5 2 - 7205 fn Q , 

486 1 -5 + 2916 ' ^ ^ 


The latter function is a new result from Ref. [24] that originates from sqq'y final states (q = 
u, d, s ). Contributions to b —» Xf'y from such final states at the NLO have been neglected in the 
previous literature because they are suppressed by phase space factors and the small Wilson 
coefficients Cz,...,g- 


Appendix D: Input parameters 

In this appendix, we collect numerical values of the parameters that matter for our branching 
ratio calculation in Section HI The photon energy cut is set to E 0 = 1.6 GeV. Our central 
values for the renormalization scales are /+ = /z c = 2.0 GeV and /j 0 = 160 GeV. 

12 Eq. (3.12) of Ref. 34] gives 4>^ A only, and contains a misprint in the coefficient at lim mc _> TOi) . 
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Masses of the b and c quarks together with the semileptonic B —> X c £u branching ratio B r p 
and several non-perturbative parameters are adopted from the very recent analysis in Ref. [491 J 1 3 I 
In that work, fits to the measured semileptonic decay spectra have been performed with op¬ 
tional inclusion of constraints from the 6-hadron spectroscopy, as well as from the quark mass 
determinations utilizing moments of R(e + e~ —> hadrons) [71]. While m c is MS-renormalized, 
m b and the non-perturbative parameters are treated in the kinetic scheme. We choose the 
option where both mb and m c are constrained by R(e + e~ —> hadrons), and m c (2 GeV) is used 
in the fit. Once the parameters are ordered as {m^kin, m c (2GeV), p 2 , p 3 D , p G , p \ s , Bdp} 
(expressed in GeV raised to appropriate powers), their central values x, uncertainties <?, and 
the correlation matrix R read [53] 


X = 

( 

4.564 

1.087 

0.470 

0.171 

0.309 

-0.135 

10.67 

a = 

( 

0.017 

0.013 

0.067 

0.039 

0.058 

0.095 

0.16 


/ 

1.000 

0.461 

-0.087 

0.114 

0.542 

-0.157 

-0.061 



0.461 

1.000 

-0.002 

-0.020 

-0.125 

0.036 

0.029 



-0.087 

-0.002 

1.000 

0.724 

-0.024 

0.049 

0.153 

R = 


0.114 

-0.020 

0.724 

1.000 

-0.101 

-0.135 

0.076 



0.542 

-0.125 

-0.024 

-0.101 

1.000 

-0.011 

-0.009 



-0.157 

0.036 

0.049 

-0.135 

-0.011 

1.000 

-0.023 


V 

-0.061 

0.029 

0.153 

0.076 

-0.009 

-0.023 

1.000 


Apart from the above parameters, the analysis of Ref. [45] serves us as a source of a numerical 
formula for the semileptonic phase-space factor 


C 


V ub 2 T[B X c eu] 
V cb Y[B ->■ X u eD] ’ 


which reads [53] 


(D.2) 


C = g(z ) {0.903 - 0.588 [a a (4.6 GeV) - 0.22] + 0.0650 [m 6ikin - 4.55] 

- 0.1080 [m c (2 GeV) - 1.05] - 0.0122 p 2 G - 0.199 p 3 D + 0.004 p 3 LS ] , (D.3) 

where g(z) = 1 — 8x; + 8^ 3 — z 4 — 12 z 2 In z and 2 : = m 2 c { 2 GeV)/m^ kin . Next, we use C in the 
expression [72j 


B S1 (E 7 > E 0 ) 


B, 


civ 


Vt*sVt b 

Vcb 


6ttem 


\P(E 0 ) + N(E 0 )], 


(D.4) 


to determine the radiative branching ratio. Known contributions to the non-perturbative cor¬ 
rection N(E 0 ) are given in terms of pp 3 Dl p 2 G and p 3 LS . The semileptonic branching ratio B c i d 
is CP- and isospin-averaged analogously to Eq. (11.3k while the isospin asymmetry effects in 
both decay rates are negligible. Thus, neither the lifetimes nor the production rates need to be 
considered among our inputs. 

13 See also the previous version [70] where more details on the method are given. 
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The remaining parameters that are necessary to determine P(Eq ) and the overall factor in 
Eq. (1D.4D are as follows: 


®em(0) 

V t *sVt b 2 

V* 


1/137.036, 

M z 

= 91.1876 GeV, 

M w = 80.385 GeV 

ra. 

0.1185 + 0.0006 

PI. 

m t , po ie = (173.21 + 0.51 + 0.71) GeV [13], 


0.9626 + 0.0012 

PI. 

— e (10,50). 
m q 


(D.5) 


For the electroweak and 0(V u b ) corrections to P(E 0 ), we also need 

ttem (Mz) = 1/128.940, sin 2 6 W = 0.23126 [13], 

M mggs = 125.7GeV p], = - 0.0080 + 0.018* [73]. (D.6) 

VtsVtb 

The quark mass ratio rrib/m q (q = u, d, s ) in Eq. (ID.51) serves as a collinear regulator wherever 
necessary. Fortunately, the dominant contributions to T(6 —)■ Xf'y) are IR-safe, while all the 
quantities requiring such a collinear regulator contribute at a sub-percent level only. They 
undergo suppression by various multiplicative factors (C 3] ... ]6 , Q d a s /7i, etc.), and by phase- 
space restrictions following from the relatively high E 0 ~ mb/3. Changing mb/m q from 10 to 
50 affects the branching ratio by around 0.7% only. We include this effect in our parametric 
uncertainty even though the dependence on m b /m q is spurious, i.e. it should cancel out once 
the non-perturbative correction calculations are upgraded to take collinear photon emission 
into account (see Refs. [M,22[Z5]J). Thus, the parametric uncertainty due to mb/m q might 
alternatively be absorbed into the overall ±5% non-perturbative error [IT] . Our range for 
m b /rn q roughly corresponds to the range [m B /mK, bn B jm^\, which is motivated by the fact 
that light hadron masses are the physical collinear regulators in our case. 

All the uncertainties except for those in Eq. (ID.ID are treated as uncorrelated. One should 
remember though that the dependence of C on a s is taken into account via Eq. (ID.3|) . 
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